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ABSTRACT 

We examine X-ray emission produced from hot gas during collisions and mergers of disk galaxies. 
To study this process, we employ simulations that incorporate cosmologically motivated disk-galaxy 
models and include the effects of radiative cooling, star formation, supernova feedback, and accreting 
supermassive black holes. We find that during a merger, the colliding gas in the disks is shock-heated 
to X-ray-emitting temperatures. The X-ray luminosity is spatially extended, rises during the initial 
stages of the merger, and peaks when the galactic centers coalesce. When a physical model for accreting 
black holes is included, the resulting feedback can drive powerful winds that contribute significantly 
to the amount and metallicity of hot gas, both of which increase the X-ray luminosity. In terms of 
their stellar kinematics and structural properties, the merger remnants in our simulations resemble 
elliptical galaxies. We find that the X-ray luminosities of the remnants with i?-band luminosities in 
the range Lb ~ 10 10 — 10 11 Lq are consistent with observations, while remnants with smaller or 
larger masses are underluminous in X-rays. Moreover, because the majority of the merger remnants 
are broadly consistent with the observed scaling relations between temperature, B-band luminosity 
and X-ray luminosity we conclude that major mergers are a viable mechanism for producing the X-ray 
halos of large, luminous elliptical galaxies. 

Subject headings: galaxies:active — galaxies:evolution — galaxies: formation — galaxies:interactions 
— mcthods:N-body simulations — X-rays: galaxies 



1. INTRODUCTION 

It is now well established that elliptical galax- 
ies contain substant i al qu antities of hot gas (sec 
iMathews fc Brighentil 120031 and references therein). 
Detailed studies of the X-ray emission produced by 
this hot gas have revealed that its luminosity is cor- 
related with the gas temperature, B-band luminos- 
ity, and stellar velocity dispersion of the host el- 
liptical gal axy llO'Sullivan. Forbes, fc Ponmai] l2001at 
IP 'Sullivan. Ponman. fc Collins! 120031) . These correla- 
tions argue for an evolutionary link between the process 
that formed the elliptical galaxy and the surrounding, 
X-ray emitting gas. 

According to the "merger hypothesis," elliptical galax- 
ies are formed when two spiral galaxies int eract and 
merge l|Toomre fc Toomrelll972HToomre1ll977D . Thus, it 
is possible that disk-galaxy interactions transform spirals 
galaxies into ellipticals and generate the ir hot gaseous ha- 
los. U sing the X-ray satellite RO S AT . lRead fc Ponmanl 
(1998) investigated the X-ray luminosit y evolution along 
the "Toomre sequence" flbom3[l927|), a series of eight 
local galaxies representing different stages of merging, in 
chronological order. Their study found that extended 
X-ray emission is produced during the encounter and 
persists as the nuclei merge. Similar results obtain 
from observations of merging and interacting systems us- 
ing the higher resolut i on ChandraX — rayObservatory 
llFabbiano et all 120011: iMcDowell et all 120031: iHuo et all 

A puzzling aspect of the analysis of iRead fc Ponmanl 
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( 1998) is that the merger remnants appear to be un- 
derluminous when compared to typical elliptical galax- 
ies, bringing into question the likelihood that merger 
remnants evolve into normal ellipticals. These authors 
suggested several methods by which a merger remnant 
could acquire an X-ray halo, late infall of tidal mate- 
rial, reacquisition of gas ejected by galactic winds during 
the merger, and ongoing mass-loss by stars born dur- 
ing a merger-induced starburst. One clue to the origin 
of the hot gas lies in the correlation between the X-ray 
and B-band luminosities. tO'Sullivan et alJ l)2001afl . also 
using ROSAT, found that Lx oc Lb for less massive 
ellipticals, consistent with their X-ray luminosity being 
primarily stellar in origin. However, the X-ray emission 
from more luminous ellipticals varies as Lx oc L B , indi- 
cating a non-stellar formation process. 

Another clue to the origin of hot gas in spheroidal sys- 
tems is its high metallicity. iHumohrev fc Buotei (2006) 
analyzed a sample of 29 X-ray luminous elliptical galax- 
ies and found that their two-phase models for the hot X- 
ray emitting plasma yielded metallicities that were typ- 
ically solar or even higher and generally correlated with 
the stellar metallicity. This result indicates that the gas 
that makes up the hot coronae surrounding ellipticals 
has been recycled from star-forming gas and likely once 
resided within the dense central regions of the galaxies. 
This also suggests that a significant amount of energy 
was injected into this gas in order to expel it from the 
depths of the potential well and redistribute it through- 
out the galactic halo. 

The work described here is a first attempt to address 
the viability of galaxy mergers in generating X-ray emit- 
ting gaseous halos consistent with observed ellipticals. 
While numerical simulations have demonstrated that 
galaxy interactions can indeed produce remnant galaxies 
whose stellar component is reminiscent of local ellipticals 
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(Barnes 1992: Barnes & Hernciuis 


?tJll99fiHHernauistll992. 


1993b: Naab&Burkert 






Robertson et all 12005: 


SDrineel et alJl2005al ICox et al.ll2005al). little is known 



theoretically about the evolution of the hot gas c ompo- 
nent in these events. iBarnes fc Hernauistl l)199fil) con- 
sidered simulations in which the gas was not allowed to 
cool radiatively and showed that in this limit it would 
be mostly shock heated to the virial temperature of the 
remnant, l eaving a co r ona o f hot, X-ray emitting mate- 
rial. Later. ICox et alJ 1)2004(1 showed that the amount of 
hot gas is dependent on the interaction orbit, i.e., more 
radial encounters produce more hot gas. Yet the ability 
of this hot gas to produce X-rays, the luminosity of the 
X-ray emission, and the correlation to the remnant stel- 
lar properties are all open questions. In particular, these 
studies did not attemp t to model all feedback processes 
in deta il. Recently, iSpringel. Di Matteo. fe Hcrnauist 
(2005bJ have developed methods for modeling feedback 
from both star formation and black hole growth in galaxy 
mergers. This addition to our galaxy modeling has signif- 
icant implications for the X-ray evolution as the feedback 
energy from an accreting black hole can generate a sub- 
stantial amount of hot gas through the production of a 
large-scale galactic wind. 

The rest of this paper is organized as follows. In § [5] 
we describe the numerical simulations upon which are re- 
sults are based. In §|3we detail the X-ray emission pro- 
duced during and after the collision of two gas-rich disk 
galaxies. We follow this by showing how the global prop- 
erties of X-rays and hot gas depend on progenitor mass 
and gas fraction in §0] § discusses our results, specif- 
ically addressing the X-ray luminosities of our merger 
remnants and the generation of hot gaseous halos ob- 
served in luminous elliptical galaxies. Finally, we con- 
clude in [Jill 

2. MERGER SIMULATIONS 

To s imulate merger s of disk galaxies, we employ GAD- 
GET2 l|Springell2005|) . which is based on a "conservat ive- 
entropy" formulation ijSpringel fc Hernauistl l2002|) of 
smoothed particle hydrodynamics (SPH) that conserves 
both energy and entropy (un like earlier versions of 
SPH; see e.g.. !Hernauistlfl993a|) . while improving shock- 
capturing. The code includes the effects of radia- 
tive cooling and heating by a uniform UV back- 
ground, and star formation, supernova feedback, and 
metal enrichment are treate d in the manner detailed in 
ISpringel fc Hernauistl (|2003) . Within this formalism, the 
sub-resolution physics of star formation and supernova 
feedback acts to maintain star-forming gas at a temper- 
ature specified by an e ffective equation of state (see e.g., 
IRobertson et alJl2004jl . As a n extension of this star for- 
mation and feedback model, ISpringel et alJ l|2005bfl in- 
troduced an additional parameter qeos, so that the ef- 
fective equation of state can be varied between that for 
an isothermal gas at T c g = 10 4 K, qeos — 0, and the 
"stiff" Springel-Hernquist equation of state, which has 
an effective temperature T c g rj 10 5 K, qeos = 1. In 
the work presented here, qeos was set to 0.25 result- 
ing in a mass-weighted temperature of star-forming gas 
~ 10 4 - 5 K. 

Our simulations also allow the presence of a supermas- 
sive black hole (BH). BHs are represented by "sink" par- 
ticles that can accrete neighboring gas at a rate given 



by the Bondi-Hoyle-Lyttleton approximation with an 
imposed upper limit equal to the Eddington rate. A 
small fraction (typically 5%) of the bolometric luminos- 
ity (O.lMc 2 , assuming an accretion efficiency of 10%) 
thermally couples to the gas surrounding the BH. These 
parameters are selected so that simulations of disk-galaxy 
mergers successfully reproduce the normalization of 
the M fjff-q relation l|Di Matteo. Springel. fc Hernauistl 
2005). The BH, when included, has a seed mass of 
w 10 5 M®, although we find that the outcome is in- 
sensitive to this choice. A mo re complete dis c ussion of 
the methodology is provided in ISpringel et alJ lj2005bf ). 

The mer gers used througho ut this paper are identical 
to those in ICox et alJ ll2005al) and ar e descr ibed in more 
detail there and in IRobertson et alJ l)2005|) . A full de- 
sc ription of the construc tion of these galaxies is given 
in ISpringel et al.l l|2005b|) : here we simply outline the 
parameter choices for the simulations discussed in this 
work. The composite system contains an exponential 
disk embedded in a dark matter halo. For simplicity we 
do not include a s pheroidal bulge. The dark-mater halo 
is initialized with a lHernauistl l(l990) profile, a concentra- 
tion c = 9, and a spin parameter A = 0.041. In Section^] 
we investigate hot gas production for mergers between 
equal mass galaxies with a range of masses, but much of 
our discussion is based on a fiducial case where the pro- 
genitor galaxies have properties similar to the Milky Way. 
The virial velocity in this case is V200 = 160 kms -1 , re- 
sulting in a halo mass of 1.4xl0 12 Mq. The exponential 
disk composes 4.1% of the total mass, with a fixed frac- 
tion / in a collisional gaseous component. Our fiducial 
simulation fixes / = 0.4. The compound galaxy is real- 
ized with 500,000 particles to represent the dark matter 
and 50,000 to represent the disk. Since the disk com- 
prises both stellar and gaseous components, a fraction 
/ x 50, 000 are consider gas and the remainder represent 
the collisionless stellar disk. 

Once the model galaxies are generated, we collide two 
disks by placing them on a parabolic orbit with an ini- 
tial separation of 143 kpc and a pericentric distance of 
7.1 kpc, resulting in a nearly radial collision. To inves- 
tigate the consequences of black hole growth and feed- 
back, we have performed simulations that include cen- 
trally accreting BHs, and are labeled as such, and others 
that do not that are referred to as "standard". Figure ^ 
shows six representative images of both types of galaxy 
mergers (i.e., with and without BHs). Details of the 
morphology, dynamics and r esulting star formation can 
be found elsewhere (see e .g. . IBarnes fc Hernauistl Il996l: 
Mih os fc Hcrnauist llflM ISpringell 1200(1 ISpringel eTaLl 
2005b; Cox et al. 200§\$). 

Figure also demonstrates the additional diffuse gas 
present when the merger includes an accreting black hole. 
Owing to the energy put into the gas by the growing BHs, 
particularly during the final merger, a significant fraction 
of the gas participates in a large-scale galactic wind. This 
"blowout" phase is very ef ficient and star forma tion is 
almost completely shut off ()Springel et alJl2005aft . leav- 
ing a rem nant that quickly red dens and evolves nearly 
passively l|Hopkins et alJl2"006a|) . Moreover, as discussed 
further in the following section, the expulsion of the gas 
during the blowout greatly enhances the X-ray luminos- 
ity relative to cases in which no BH is present because 
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Fig. 1. — Time sequence of the projected gas density, shown in blue scale, with X-ray contours overlaid (in black). The top six panels 
demonstrate the collision for our standard, no-black hole merger, while the bottom six panels are the case in which a central supermassive 
black hole is included. Each panel is 100 kpc square and the color and X-ray contour scales are identical. The simulation time, in Gyr, is 
indicated in the top left of each plot, and the merger occurs at T fs 1.6 Gyr. 
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not only is there additional gas throughout the halo, but 
this gas is highly metal-enriched. 

While a complete survey of the parameter space that 
describes the merger of two disk galaxies is beyond the 
scope of this work, it is important to understand how 
the production of hot gas depends on our initial assump- 
tions. To this end, we merge disk galaxies with varying 
gas fractions and masses and describe the outcomes in 
Section H 

3. A CASE STUDY 

In this section we provide a detailed look at the evo- 
lution of hot gas and X-rays during one merger. While 
the X-ray luminosity, gas temperature and mass depend 
strongly on the initial halo size/mass and disk gas frac- 
tion (discussed in §UJ|, the production of hot gas during a 
galaxy collision is a general feature of all our simulations. 

3.1. X-ray emission 

Our analysis of the simulations assumes that X-rays 
are produced by the cooling of hot, diffuse gas. In real- 
ity, the observed X-ray emission has contributions from 
stellar sources such as supernova remnants and X-ray 
binaries, as well as accreting black holes and the hot 
phase of the interstellar medium. Because the compos- 
ite emission from these sources is generated mostly in 
the dense central regions of the galaxies, this luminosity 
may be heavily extinguished owing to the large column 
de nsity of intervening gas and du st. In fact, as described 
bv lHo'pkins elaD l)2005dt l2006a|) . we find that the cen- 
tral accreting BHs in merger simulations remain obscured 
for the majority of their active lifetimes. To simplify our 
analysis, we restrict ourselves to X-ray emission from the 
diffuse hot gas since it should be relatively unaffected by 
obscuration. 

We follow [Navarro et alJ l)1995j) and assume the X-ray 
luminosity for each SPH particle can be estimated from 
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where m p is the proton mass, /i is the mean molecular 
weight (0.6 for a fully ionized primordial plasma), and 
m g&s,i, Pi, and Tj are the mass, density and temperature 
of the i th gas particle in cgs units, respectively. Equa- 
tion JQ| assumes that the primary mechanism for X-ray 
emission is thermal bremsstrahlung, an assumption con- 
sistent with the zero-metallicity cooling included in the 
simulation. However, metal enriched gas with a tem- 
perature of ~ 10 6 K cools primarily through metal-line 
emission, a mechanism that is more efficient and would 
thus produce a much larger X-ray luminosity. In this 
sense, the X-ray emission computed via Equation Q is 
a lower limit. 

In light of the increased X-ray luminosity expected 
from metal-enriched gas, we supplan t equation l|T|) with 
X-ray emission calculated using a IRavmond fe Smith 
l|1977l hereafter RS77) code that does include metal-line 
emission. The gas metallicity is provided by t he star- 
formation model fsee lSpringel fc Hern auist 2003) assum- 
ing a metal yield of 2%, and we include emission from 
0.1 to 2 keV. 

The bolometric luminosity is then calculated according 
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Fig. 2. — Bolometric X-ray luminosity owing to diffuse 
hot gas during simulated major mergers. Shown are mergers 
with ("BH", blue) and without ("std", red) accreting black 
holes. Each simulation has two curves, one, in which the 
gas is assumed to have zero metallicity (blue, solid line, and 
red, asteri sks), and, the other, in which the emission is calcu- 
lated using Ravmond & Smith ( 1977, RS77) models, which include 
metal-line cooling (blue, open squares, and red, filled circles). 
The arrow at T = 1.5 Gyr shows the time of the final merger 
in both simulations. 
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where the sum is carried out over all gas particles with 
density below the critical density for star formation 
Pthresh ~ 10~ 2 Af Q pc and a temperature greater than 
10 5 ' 2 K. Gas particles that meet these criteria are fully 
ionized (for zero metallicity) and are deemed "hot" or 
"X-ray" gas from hereafter. We note that our low tem- 
perature (10 52 K) threshold for gas to be included in 
the X-ray luminosity probes temperatures below current 
observational limits. However, this cut-off does serve as 
an accurate gauge of the amount of hot, ionized gas. 

3.2. Time evolution 

Figure |3 shows the bolometric X-ray luminosity as de- 
termined by equation J2J for the entire major merger. 
For both simulations, two X-ray luminosities are shown. 
One, where the X-ray luminosity is calculated via equa- 
tion , in essence when the gas is assumed to have zero 
metallicity, and the second uses the RS77 code to include 
metal-line cooling. As expected, the inclusion of metal- 
licity dependent cooling greatly enhances the X-ray emis- 
sion produced during the galaxy merger. In addition, the 
spatial distribution of the X-ray emission is shown with 
overlaid contours in Figure ^ 

Initially, the X-ray emission from both galaxies is neg- 
ligible because we ignore contributions from the star- 
forming interstellar medium. X-ray emission begins 
when the gas disks first interact, in shocks that lie di- 
rectly between the two disks (see T = 0.28 Gyr panel 
in Fig. [TJ. In general, the majority of the X-ray emis- 
sion is produced by shock-heated gas, and thus the X-ray 
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emission essentially tracks shocks that occur during the 
galaxy merger. 

Figure n also shows that when the galaxies separate, 
the X-ray emission is concentrated around the two galac- 
tic nuclei. This emission resembles that in observed 
systems t hat are in a similar dynamical state, such as 
the Mice l)Readll2003^ . Figure ^ also demonstrates that 
when BHs are present, additional hot gas is generated 
by thermal feedback from BH accretion beginning at 
T = 0.7 Gyr. This results in an Lx that is over an 
order of magnitude larger than the no-BH case. 

The excess X-ray emission from T = 0.7 to 1.5 is gen- 
erated by small-scale winds in the still-separated disks. 
These winds are, in general, collimated and escape per- 
cndicular to the disk plane, similar to those shown in 
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ringel et alJ l)2005bl their Fig. 9). In the simulations 
presented here, the winds are generated by the energy in- 
put by the accreting black holes and thus are present only 
in the BH simulations. Collimated winds such as these 
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One of the most prominent features in Figure El is the 
coincidence of the peak X-ray luminosity with the final 
merger, regardless of whether or not BHs are present. 
This X-ray emission results from the abundant shocks 
that occur as the gas disks attempt to follow the collision- 
less stellar nuclei. Gas is heated to 10 6 — 10 8 K and forced 
out of the central regions owing to shocks and pressure 
forces. Hence, the X-ray luminosity is spatially extended, 
as can be seen in Figure ^from T ~ 1.5 — 1.9 Gyr. 

The maximum X-ray emission ranges from 8 x 10 39 to 
2.5 x 10 41 ergs s _1 , depending on the presence of a BH 
and if metallicity-dependent cooling is included. These 
luminosities are consistent with observations of the X- 
ray emission seen in i nteracting and merging systems 
(iFabbiano et all 120011: iXia et all 120021 iMcDowell et al l 
1 20031: IHuo etaD 2001 iJenkins et all 120051: iReadl 120051: 
ISmith et alf2 005) as well as tha t expected on theoretical 
grounds ll,Tog fc Solomon1ll992^ . Because the X-ray lu- 
minosity is a direct consequence of the collision between 
gas in each progenitor disk, we expect a range of peak 
luminosities depending on the gas content in each galaxy, 
their (relative and cumulative) mass and the merger or- 
bit. The mass and gas fraction dependence of Lx is 
investigated further in § 0] 

After the galactic nuclei coalesce, at T ~ 1.6 Gyr, 
the X-ray luminosity decreases. When no BH is present, 
Lx drops slowly as gas settles into hydrostatic equilib- 
rium, and tidal material slowly returns to the remnant. 
However, when an accreting BH is included, the drop in 
Lx is more pronounced owing to the "blowout" phase 
which expels a significant amount of metal-enriched, X- 
ray emitting gas. The eventual fate of gas caught up in 
the galactic wind is dictated by its energy relative to the 
depth of the galaxy potential well. In the simulation we 
describe in this section, a substantial amount of gas is in- 
jected with enough energy to completely unbind it from 
the remnant. However, even if a significant amount of gas 
is ejected from the system, the blowout does not leave the 
remnant completely devoid of gas. A substantial amount 
(~ 75%) of the remaining gas remains bound to the rem- 
nant galaxy, is spread throughout the dark matter halo, 
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Fig. 3. — Time evolution of gas mass during the disk galaxy 
merger. Shown with a solid black line is the total gas mass. Dur- 
ing the simulation, gas is steadily depleted owing to efficient star 
formation. The (blue) line with crosses is the mass of gas which 
emits X-rays. The (red) line with asterisks is mass of gas which 
becomes unbound from the system. Two plots are shown, one with 
and the other without a centrally accreting black hole. 



and makes up the X-ray corona. 

Figure [21 provides an accounting of the gas mass for 
mergers with and without accreting BHs. This figure also 
shows the mass of unbound (i.e., E > 0) gas and the mass 
of gas that contributes to the X-ray luminosity. Gas that 
is not unbound or X-ray emitting is typically cold and, 
depending on its density, star formin g. Consistent with 
the findings of ISprineel et alJ l|2005a|) there is little cold 
gas in the remnant that includes a BH and star formation 
is efficiently quenched subsequent to the merger at T ~ 
1.6 Gyr. 

Figure |3 demonstrates that the BH has a significant 
impact on the amount of X-ray and wind gas. Owing to 
the merger-induced energy input from the BH, and the 
significant wind that ensues, the BH system has over 10 
times the mass of unbound gas as compared to the no 
BH case. The BH feedback also contributes the X-ray 
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gas and results in nearly 3 times the amount of mass 
as compared to when no BH is present. At first the 
differing amounts of X-ray gas in the two cases seem at 
odds with the identical X-ray luminosity (for emission 
by thermal bremsstrahlung) demonstrated in Figure [21 
However, as shown in § 13.51 (see Fig. |HJ) the gas within 
the inner 10 kpc of the remnant without a BH is over 
an order of magnitude denser than the no BH remnant 
and thus is more efficient at producing X-rays. We note 
that the thermal structure is quite similar between the 
two remnants, as is shown in § 13.41 

The hot gas in the remnant relaxes on a timescale that 
is short when compared to that of the stellar component 
of the halo, as reflected by both the regularity of the 
X-ray isophotes (see the bottom-right panel in Fig. ^| 
and the steady X-ray luminosity after T = 3.0 Gyr, 
a little more than 1.5 Gyr after the merger. In con- 
trast, stellar material in tidal tails is accreted by the 
remnant to pro duce long-lived shells, loop s, and other 
fine struct ure iHerna uist fc Snergell 1 19921) bv "phase - 
wrapping" l)Quinnll984HHernauist & Quinnll988lll989D . 
Similar features are seen around many ellipticals 
HTvTalm fc Carter! IT983 IsThw^el 1798(1 ISchweizer et al.l 



1990). Because this debris is on nearly radial orbits, gas 
in the tails is stripped out an d deposited in the halo 
or inner regions of the r emnant ijHernauist fc Weillll992t 
iWeil fc Hernauisrt lT993?) . and the shells quickly redden. 

The BH-driven wind is also efficient at dispersing met- 
als. Because the blowout is triggered by the centrally 
located BH, gas that is ejected from the system has been 
significantly enriched by metals. To demonstrate this 
more explicitly we refer to Figure 0] Plotted are con- 
tours of the gas density at T w 0.9 Gyr, shortly after the 
first passage of the two disks. Overlaid are {red) dots 
that denote the gas particles that will eventually con- 
stitute the hot X-ray emitting gaseous halo. When no 
BH is present, the gas in the X-ray halo was once almost 
exclusively composed of loosely bound tidal material. In 
these low-density regions, little star formation occurs and 
thus metal enrichment is only modest. While this same 
loosely bound tidal material will make it into the X-ray 
halo of the BH remnant, the majority of the X-ray gas 
in this case comes from gas that was once more centrally 
concentrated. This gas is spread throughout the disk and 
has been significantly enriched by star formation. 

For the BH merger shown in Figure X-ray gas in the 
remnant has an X-ray weighted mean metallicity of 0.9 
solar. The unbound wind material has an even higher 
metallicity, nearly twice solar. In this sense, the wind 
could be considered "mass loaded." The high metallic- 
ity of the gas suggests that the wind is made up pri- 
marily of the densest, most metal-rich gas. In contrast, 
without the energy input from a BH, the X-ray gas is 
1/10 Zq, and the small amount of unbound gas is less 
than 1/20 Zq. Without the energy input from the BH 
essentially the entire metal budget is locked in stars or 
star-forming gas. 

The dispersal of metals by the wind also has a profound 
influence on the X-ray emission. Because the metallicity 
of the X-ray-emitting gas is nearly 10 times larger in 
the BH remnant as compared to the standard remnant, 
the X-ray luminosity is greatly increased. We show the 
X-ray luminosity when metal-line cooling is considered 
in Figure |3 We hereafter restrict our analysis to X-ray 
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Fig. 4. — Contours show the gas density subsequent to the first 
passage of the interacting dark galaxies. Tidal tails and a bridge are 
clearly present. Overlaid are (red) dots that denote gas particles 
that eventually end up as part of the hot, X-ray emitting gaseous 
halo. Two images are shown, one with and the other without a 
centrally accreting black hole. 



emission calculated using the metal-line cooling of RS77. 

3.3. Time evolution of Lx/Lb 

Because galaxies encompass a range of masses, studies 
of the X-ray emission caused by galaxy interactions have 
typically normalized the X-ray emission by the .B-band 
luminosity of the host galaxy ijRead k, Ponmanl Il998t 
IQ'Sullivan et al.M2001b|) . In order to facilitate a more 
direct comparison to these studies, we show Lx/Lb for 
our mergers in Figure [5j For simplicity, we ignore the 
effects of extinction when computing L b ■ This approxi- 
mation is likely to be valid during the post-merger phases 
but may dramatically overestimate the luminosity dur- 
ing the highly obscured merger stage. In any case, many 
of the features in Lx/Lb are very similar to those of 
just Lx (shown in Fig. |2J). This result is not surprising 
given the moderate gas fraction of our progenitor disks, 
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Fig. 5. — As in Fig. |21 except the X-ray luminosity has been 
divided by the _B-band luminosity. For clarity, only the metal-line 
X-ray emission is plotted. Again, the arrow at T = 1.5 Gyr shows 
the time of the final mer ger in both simu lations, and the slope 
shows the observed slope IQ'Sullivan et alJI2001bl. OFP01). 



and hence the modest fraction of stellar mass that is pro- 
duced during the merger. For the most part, Lb slowly 
decreases over time owing to an aging stellar population. 

Both Lx and Lx/Lb (Figs. |21 and 03 respectively) 
demonstrate a characteristic rise and fall of X-ray emis- 
sion. The peak of this emission is coincident with 
the merger, a trend t hat is also seen observationally 
ijRead fc Ponmanll 19981 their Fig. 13). One discrepancy 
may be the decline from the peak which is larger in the 
simulations than the observations. However, Figure El 
shows a wide range of peak luminosities depending on 
the presence of a BH and the metallicity of the gas. The 
luminosities will also fluctuate owing to the gas content of 
the progenitor disks as well as the merger orbit, compli- 
cating a quantitative comparison with the observations. 

We can also compare the X-ray luminosity from the 
simulations subsequent to the merger to the s tudy of 
disturbed ellipticals bv 10 'Sullivan et all <|2001b|) . These 
authors determined that L x / L B steadily increases with 
remnant age when age is measured by th e fine structure 
parameter £ l|Schweizer fc Seitzerlll99^) and is zero at 
the last merger event. The slope of the best-fit relation 
by O' Sullivan et al. is shown as a short line in Fig- 
ure ("OFPOl"). The shallowness of this X-ray evolu- 
tion demonstrates that the increase of Lx I Lb over time 
is a very subtle effect and that the hot gaseous halos of 
elliptical galaxies are grown over long periods of time. 

For the simulation that includes BHs, the Lx/Lb ra- 
tio is nearly constant after the merger, a by-product of 
the equivalent fading of both Lx and Lb- While we 
followed each remnant for only a short period of time af- 
ter the merger and the isolated nature of the galaxies in 
our simulations make the long term luminosity evolution 
fairly uncertain, it does not appear that our remnant 
Lx/Lb increases at the same rate as observed merger 
remnants. 

3.4. Remnant emission 
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Fig. 6. — Azimuthally averaged X-ray surface brightness profile. 
The line types are identical to Fig. [F] 



As discussed in §^the "merger hypothesis" posits that 
the interaction and merger of two spiral galaxies leads 
to the formation of an elliptical. In general, our simula- 
tions support this picture as we are left with a spheroidal 
stellar component surrounded by a halo of hot, X-ray- 
emitting gas, similar to observed elliptical galaxies. In 
this section we investigate this in more detail by directly 
comparing the X-ray emission of our remnants to obser- 
vations of post-merger and relaxed elliptical galaxies. 

To begin, we show the X-ray surface brightness in Fig- 
ure The emission is calculated using the RS77 metal- 
line cooling and is azimuthally averaged from 100 random 
viewing directions. The presence of BHs has a significant 
impact on the shape and radial extent of the X-ray emis- 
sion. Apparently, energy injected by the BH decreases 
the emission on small scales (< 3 kpc), and increases 
the emission on large scales (> 3 kpc). In both cases, 
however, there is detectable emission to ~ 100 kpc. 

It is common practice to fit the observed X-ray surface 
brightness profiles of elliptical galaxies to a beta model 
of the form 

EX (R) = SX0 [1 + (R/Rcore) 2 } ^ +0 ' 5 , (3) 

where Ex.o is the central surface brightness, R cole is the 
core radius, and /? quantifies the slope at large projected 
radii. We fit equation © to the surface brightness of 
our remnants and find the best fit to yield (3 — 1.85 and 
-Rcore = 43.3 kpc for the BH remnant and (3 = 2.12 and 
-Rcore = 38.6 kpc for the remnant without a BH. 

We note that the above fit was performed over the en- 
tire surface brightness profile, which includes emission 
out to ~ 100 kpc. At these radii, the X-ray surface 
brightness is over 5 orders of magnitude fainter than the 
central regions. This large dynamic range is nearly 3 or- 
ders of magnitude more than is probed by the standard 
integration time of observations. Thus, a more realistic 
comparison can be achieved by restricting our fits to radii 
where the X-ray surface brightness is within 2 orders of 
magnitude of the peak. This corresponds to 12 and 50 
kpc for the standard and BH remnant, respectively. In 
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Fig. 7. — Projected temperature and metallicity profiles for 
the merger remnants. Both temperature and metallicity are X-ray 
luminosity weighted and azimuthally averaged. Line types are as 
in Fig. [5] 



this case, we find f3 = 0.80 and i? C ore = 18.4 kpc for the 
BH remnant and /3 = 0.53 and R cmo = 0.73 kpc for the 
remnant without a BH. 

The /3-profile fits restricted to the inner surface bright- 
ness profiles are in much better agreement with obser- 
vations than the unrestricted fits. Observed merger 
remna nts have m e asure d values between 5 and 
1.15 (INolan et a J l200l IFabbiano fc Schweizerl H995t 
iKim fe Fabbianol 12003) . We also note the much flatter 
profile and larger core radius of the BH remnant. These 
features are observed in sever al X-ray faint early-type 
galaxies l)Q 'Sullivan fc PonmanlfeOO^ . 

In addition to the X-ray surface brightness profile, an 
increasing number of galaxies have well determined tem- 
perature and metallicity profiles. Acquiring this infor- 
mation has been possible with the high spatial and spec- 
tral resolution of Chandra and XMM — Newton, as 
well as the proper accounting of point sources and com- 
plex modeling of the hot X-ray emitting plasma. In par- 
ticular, recent work has shown that multi-temperature 
plasma models provide a much better fit to the data 
and result in a larger e stimate for the hot gas metallicity 
l|Buote fc Fabianlll998T) . 

Figure shows the temperature (top) and metallicity 
(bottom) profiles for our merger remnants. Regardless of 
the presence of BHs, the temperature is approximately 
isothermal within ~10 kpc, and a strongly decreasing 
function at larger radii. The additional energy injected 
by the BH uniformly heats the gas at all radii without 
affecting the slope. However, inside of ~ 2.5 kpc, the 
temperature difference is much larger owing to different 
cooling times of the hot gas. Without a BH, gas has a 
shorter cooling time and a cool(er) core has formed. In 
contrast, the energy input by the BH generates a hot(ter) 



cusp. We note that these curves are expected to have 
some time variability since gas has cooling times less than 
the Hubble time in both instances (see § I3.5fl . 

In comparison to the temperature, the differences in 
the metallicity profiles are much more striking. As de- 
scribed earlier, the BH induced galactic wind is "mass- 
loaded" and transports a significant amount of metals to 
large radii. This metal-rich wind results in an increasing 
metallicity profile. The largest radial bin is supersolar 
and the most metal rich. As we showed in § 13.21 the 
trend of increasing metallicity with radius continues to 
the unbound wind material as well. Without a BH, met- 
als are retained very near to where they are produced and 
hence the metallicity profile is flat or slightly decreases. 

The majority of the spatially resolved data exists for 
X-ray luminous early- type galaxies. These systems typ- 
ically have nearly isothermal temperature profiles and 
metallicity gradients, although there is very often a 
slight de crease in temperature at small or large radii, 
or b oth llO'Sullivan fc Ponmanl l2004alU: iSivakoff et alJ 
l2004UKim fc Fabbianoll2003HMushotzkv et al.ll994j) . At 
first glance, the profiles presented in Figure appear to 
be at odds with the observed trends. However, as we 
mentioned previously, restricting our analysis to the in- 
ner regions would render the temperature profile much 
closer to isothermal. The metallicity gradient may be 
indicative of mixing between metal-rich gas produced 
during the merger and diffuse low-mctallicity gas that 
is accreted by the galaxy, a process that is not included 
in the present simulations. We also note that X-ray lu- 
minous galaxies arc generally large ellipticals, and thus 
our comparison is best directed toward merger remnants, 
which, because they are X-ray faint, do not usually have 
measured temperature or metallicity profiles. 

One generic feature of observed ellipticals and 
merger re mnants is the nea r solar metallicity of X- 
rav gas iNolan et alJ 120041: IKim fc Fabbianol l2003t 
iHumphrev fc Buotd 12006^ This is generally repro- 
duced in the simulations that include a BH. Because 
the no-BH merger produces an X-ray halo of low metal- 
licity, this suggests a mechanism to produce the low 
metallicities observed in a f ew X-ray faint ellipticals 
(|Q'Sullivan fc Ponmanl2004bJ) . This could occur if black 
hole growth is inefficient, as would be the case if the 
progenitors are gas-poor, or if the emission is so faint 
that one could only detect the central part of the rising 
metallicity profile present in the BH remnant. 

3.5. Spherical Profiles 

Figure |S] shows the spherically averaged gas density 
(top), cooling time (middle), and entropy (bottom) for 
X-ray gas and provides a general picture of the state 
of hot gas in the merger remnants. While the remnant 
profiles are quite similar regardless of the presence of 
BHs, there is a distinct transition at ~ 10 kpc. Inside 
of this radius, gas in the no-BH remnant has a higher 
density, and thus a significantly shorter cooling time, and 
lower entropy. Outside of this radius, the trend reverses, 
mirroring what was seen in the X-ray surface brightness 
profiles in Figure |SJ 

It is noteworthy that the temperature and metallic- 
ity profiles shown in Figure do not show a break. In 
this case, the BH remnant has a higher temperature and 
metallicity at all radii. Together, these trends suggest 
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Fig. 8. — Spherically averaged density, cooling time, and entropy 
for X-ray emitting gas. Line types are as in Fig. [S] 



that feedback from the BH has pumped energy and met- 
als into the gaseous component and affected diffuse gas 
at all radii. In response to this injection of energy, gas 
flows outward (and some even escapes in a wind), fill- 
ing more of the dark matter halo than when no BH is 
present. 

Even though gas that is located at > 10 kpc is rela- 
tively similar in temperature and density in both the BH 
and no-BH cases, there is a significant difference in the 
cooling times, caused by the higher metallicity of gas in 
the BH remnant. In fact, much of the existing X-ray halo 
would cool within the next several Gyr in this model. 

Finally, the spherical gas profiles display kinks at ~ 
3 kpc. While very pronounced in the cooling time and 
entropy profiles, this kink exists in some form in all the 
profiles, both projected and spherical. Features such as 
this are quite common in the simulations with BHs and 
result from a small number of residual cold gas clumps 
within the inner few kpc. 



4. SCALING WITH PROGENITOR GAS FRACTION AND 

SIZE 

Because hot gas is generated by shocks that attend the 
disk-galaxy collision, and is bolstered by feedback from 
black hole accretion, it is natural to expect that the X- 
ray emission will depend on the amount of gas present in 
the progenitor disk galaxies as well as their masses. To 
understand this relationship better, we perform a series 
of merger simulations that span a range of progenitor gas 
fractions and masses as well as the orbits by which the 
disk galaxies merge. 

All the progenitor disk galaxies used for this parameter 
exploration have identical halo properties such as concen- 
tration, spin parameter, and mass fraction. While hold- 
ing these properties constant, we systematically change 
the disk gas fraction or virial velocity, or both. In 
total, we construct three halos smaller (V200 = 56, 
80, and 115 km s _1 ) than our fiducial model (V200 = 
160 km s _1 ) and three that are larger (V200 = 225, 320, 
and 500 km s _1 ). These virial velocities correspond to 
total masses spanning 5.8xl0 10 - 4.2 xlO 13 M@, a range 
of nearly 3 orders of magnitude. For each galaxy the ex- 
ponential disk composes 4.1% of the total mass, with a 
fixed fraction / in a collisional gaseous component. The 
compound galaxies are realized with 500,000 particles to 
represent the dark matter and 50,000 to represent the 
disk. Since the disk comprises both stellar and gaseous 
components, a fraction / x 50, 000 are considered gas and 
the remainder represent the collisionless stellar disk. In 
what follows we investigate four values of /, 0.05, 0.2, 
0.4, and 0.8. Black holes are included in all of the simu- 
lations described in this section and the X-ray emission is 
calculated using the metallicity-dependent RS77 model. 

Before we begin our exploration of the effects of varying 
the progenitor mass and gas fraction, we wish to deter- 
mine if the fully co-planar progenitor disk orientations 
and radial orbit of our fiducial disk-galaxy merger signif- 
icantly affect the resultant X-ray luminosity, discussed at 
length in § |3 To this end, we merged the identical pro- 
genitor disks on the fiducial orbit with three alternate 
disk orientations. Besides our fiducial co-planar orien- 
tation (0i = 2 = 0i = 0-2 = 0), we also simulated a 
tilted orientation (0i = 30, 81 = 60,02 = —30, #2 = 45), 
another tilted orientation (0i = — 109,6*1 = —30,02 = 
71,02 = —30), and a polar orientation (0i = 60, &i = 
60,02 = 150,6*2 = 0), where and 9 are the standard 
spherical coordinates, and the subscripts 1 and 2 denote 
the two progenitor disks. The X-ray evolution for these 
runs is presented in the left panel of Figure El In the 
right panel of this figure, we show the X-ray evolution 
for our fiducial progenitor disk and fiducial orientation 
merging on six orbits with different angular momenta, 
one with less and five with more. 

Together, the orbits and orientations presented in Fig- 
ure El display a range of X-ray luminosities. At any par- 
ticular time, Lx may differ by over an order of magni- 
tude. Much of this variability is attributable to the effect 
that orbits and orientations have on the consumption of 
gas, and thus the effective gas fraction at the time the 
galaxies interact. An additional factor is the growth of 
the black hole, as larger black holes will input more en- 
ergy and produce additional hot gas. 

To understand how the progenitor gas fraction influ- 
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Fig. 9. — X-ray luminosity calculated using the mctallicity-dcpcndcnt RS77 model for several disk galaxy major mergers. In all cases, 
the progenitor disks are identical, while their orientation with respect to the orbital plane (left), or their merging orbit (right) has been 
varied. 
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Fig. 10. — X-ray luminosity (RS77) during a major merger 
between disks with varying progenitor gas fractions. The color 
and thickness of each line is proportional to the gas fraction, with 
lighter, thicker lines representing the high gas fraction mergers and 
dark thin lines the low gas fraction mergers. 



Fig. 11. — Shown is the X-ray luminosity (RS77) during 
disk galaxy major mergers for progenitors of varying mass, from 
5.8xl0 10 Mo to 4.2xl0 13 M Q . The mergers between more mas- 
sive galaxies produce larger X-ray luminosities and these curves are 
drawn with progressively lighter, and thicker, lines. 



ences the X-ray luminosity without the complications as- 
sociated with various orbits and orientations, we resimu- 
late the fiducial merger while systematically varying the 
progenitor gas fraction. The resultant X-ray evolution 
is presented in Figure ^| As the progenitor gas frac- 
tion increases, the peak X-ray luminosity also increases. 
There is also an increase in Lx prior to the merger, from 
T = 1.0 to 1.8 Gyr, owing to the larger BHs and addi- 
tional hot gas produced by their feedback. Because the 
BH feedback drives a powerful wind in the high gas frac- 
tion mergers, the X-ray luminosity drops subsequent to 
the peak. There is very little difference in the remnant 



X-ray emission for gas fractions greater than 20%. 

For the lowest gas fraction progenitor in Figure fTOl the 
X-ray emission is very small prior to the final merger. 
Many of these trends with gas fraction can be understood 
in terms of the feedback provided by the BH driving a 
large scale wind. As the gas fraction (and black hole 
mass) increase the wind becomes dominant and affects a 
larger amount of gas. Thus, the temperature and entropy 
of the remnant gas steadily increase in remnants with 
larger progenitor gas fractions. These effects will be the 
subject of future work that details the winds produced 
by black hole feedback (Hopkins et al., in preparation). 



X-Rays from Hot Gas in Mergers 



11 



The last dependency we explicitly address is progeni- 
tor mass. For this, we merge galaxies with a variety of 
masses on orbits that are self-similar, i.e. the pericen- 
tric distance is scaled in proportion to the size of each 
disk. Figure^Jpresents the X-ray evolution during these 
mergers. In all instances, the progenitor gas fraction was 
80%. 

There is a regular progression for more massive pro- 
genitors to produce more X-ray emission. This is to be 
expected, if gas settles into hydrostatic equilibrium in 
each remnant halo. There is also a subtle difference be- 
tween the maximum X-ray luminosity and the remnant 
X-ray luminosity. Both the most and least massive galax- 
ies have relatively little difference between the maximum 
and remnant Lx , while moderate masses have a large dif- 
ference. This trend results from two competing effects. 
First, the impact of black hole fe edback varies with sys - 
tem mass in a differential manner ijHopkins et al.l20 06b'). 
being increasingly more violent for more massive galax- 
ies. Second, the gas consum ption is much more effi cient 
for more massive systems ( Roberts on et alJ l2005j) and 
thus they have lower effective gas fractions. 

Finally, we summarize the production of X-rays and 
hot gas during the various disk-galaxy major mergers by 
compiling a number of global quantities for each merger 
remnant. In order to characterize the X-ray luminosity 
we determine the peak and remnant Lx for each galaxy 
merger. We also find the mass of gas that composes 
the hot X-ray halo and the mass of gas that becomes 
unbound from the galaxy potential owing to the black 
hole feedback-induced blowout. These global quantities 
are shown in Figure lT2l for our entire merger series, which 
includes the seven different mass initial conditions, four 
initial disk gas fractions, and four disk orientations. The 
scatter resulting from the four initial disk orientations 
simulated is shown by the la error bars for all quantities. 

The top panels of Figure El demonstrate what was 
also apparent in Figure 1111 namely, that more massive 
galaxy mergers produce more X-ray luminosity. For a 
self-similar model we expect Lx ~ c 4 , which is close 
to what is present in the simulated remnants. There is, 
however, a slight flattening at high a, which reflects the 
more efficient star formation and thus less relative gas in 
these systems. 

Apparent in the middle panels of Figure El is a trend 
for low-mass mergers to produce more unbound material. 
This is opposite to the trend found for X-ray gas mass, 
where the most massive systems have a larger fraction. 
Both trends hold across all progenitor gas fractions in- 
vestigated. However, there is a systematic increase in 
the fraction of unbound gas with progenitor gas fraction. 
This trend reflects the increased wind efficiency induced 
by the black hole feedback for these systems. 

The bottom panels of FigurclT2lshow the metallicity of 
the X-ray and unbound gas. In nearly all the remnants, 
both the unbound and X-ray gas have solar or superso- 
lar metallicity. There is a slight decrease in metallicity 
at high masses, as th ese systems have very s hort gas con- 
sumption timescales l)Robertson et al.l2005j) . Thus many 
of the metals become locked up in the stars, producing a 
correlation between remnant mass and stellar metallicity. 

5. DISCUSSION 



The previous section demonstrates that galaxy mergers 
produce significant X-ray emission and that each merger 
remnant contains hot, X-ray emitting gas. However, 
we have not addressed whether this X-ray emission is 
comparable to observed ellipticals. Observations indi- 
cate that merger remnants are underluminous when com- 
pared with elliptical galaxies and that ellipticals build 
their X-ray emission slowly, over the course of several 
Gyr ijO'Sullivan et alJl2001bft . In this section we com- 
pare our merger remnants to common scaling relations 
for elliptical galaxies and speculate on their future evo- 
lution. 

5.1. The Lx-Lb relation 

One of the most common relations used in the study 
of X-ray emission from individual galaxies is the Lx- 
Lb relation. Figure fHfl shows our entire series of gas 
fraction and size/mass remnants (from Fig. I12f> in the 
Lx-Lb pl ane. Also shown on th is plot as points are 
data from IQ'Sullivan et al.l i2001atl and luminous ellip- 
ticals from 10 'Sullivan et alJ l)2003|) as crosses. The dot- 
dashed line (with the shallower slope) that demarcates 
the lower bound of observed galaxies is an estimate of 
the po int source contrib ution to the total X-ray lumi- 
nosity (jCiotti et al.lll99l|) . Galaxies that reside near this 
line have X-ray emission that is presumably dominated 
by point sources. Shown with a dashed line of slope 
~ 2.7 is the best-fit rela tion for luminous ellipticals from 
IQ'Sullivan et~aTl ll2003h . 

Our merger remnants trace out an arc that passes 
through the observed galaxies. Because we have not in- 
cluded any contribution to the X-ray luminosity from 
point sources, the low-mass remnants fall well below the 
observations. This is consistent with the expectation 
that these systems radiate little diffuse gas emission or 
that the diffuse gas is too cold to be easily detected. Our 
moderate Lb (~ 10 10 — 10 11 ) galaxies are very compara- 
ble to observed ellipticals, while the higher Lb systems 
are underluminous. 

The comparison between observed ellipticals and our 
remnants is muddled by two facts, one is the uncertain 
L b owing to the unknown star formation history of stars 
that are formed before we begin our simulation, and the 
second is the baryon fraction of our initial disk galaxies, 
which we assume is less than one-third of the cosmic 
mean. 

To address the first point, we assume that stars present 
at the beginning of the simulation have been steadily 
formed during the 5 Gyr prior to our T = 0. Metallicities 
are selected randomly between solar and 1/10 solar and 
is calculated using the models of iBruzual fc Oharlotl 
(2003). We also note that the stopping point for our 
simulations is arbitrary. An additional 5 Gyr of evolution 
would decrease Lb by another factor of 2. 

The second point above, the initial baryon fraction, 
is particularly relevant for our most massive systems, 
whose closest analogs are probably groups rather than 
individual galaxies. Groups tend to h ave baryon frac- 
tions closer to the cosmic mean value ({Sanderson et alJ 
2003) and thus it may not be surprising that these rem- 
nants are underluminous. We also find it intriguing that 
X-ray emission is only dete cted in groups that h ave at 
least one early- type galaxy ijMulchaev et al.l l2003). This 
correlation naturally arises when the X-ray emission and 
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Fig. 12. — Global properties of the hot gas produced during disk-galaxy major mergers which include an accreting BH and its associated 
feedback. Each column shows results for mergers with identical progenitor gas fractions, indicated in the top panel of each column, for a 
variety of progenitor masses, as measured by the central stellar velocity dispersion of the remnants. The top row shows the maximum and 
remnant (RS77) X-ray luminosity during the merger. The middle row shows the fraction of the progenitor gas that is either physically 
unbound (i.e., E > 0) from the remnant or that constitutes the hot X-ray halo. The bottom row shows the average metallicity of the 
unbound wind and X-ray gas. 



early-type galaxy are formed in a common process, as we 
suggest here. 

5.2. The Lx-Tx relation 

Figure 1141 shows the remnant X-ray luminosity ver- 
sus X-ray-weighted average temperature for our merger 
remnants in the Lx-Tx plane. The Lx-Tx relation is 
commonly used to study X-ray emission from clusters, 
groups, and more recently, elliptical galaxies. The best- 
fit relation for ellipticals appears to have an offset in 
both slope and normalization from the relations found 
for groups and clusters (which themselves are not iden- 
tical), indicating that these systems have had an injec- 
tion of energy (or entropy) above that exp ected from 
purely self-similar gravitational collapse fsee iVoit] 120051 
for a review). On this figure we also plot the data of 
IO 'Sullivan et alj <|2003f) . which consists primarily of lu- 
minous ellipticals, as crosses, and their best-fit relation 
which has a slope of 4.8. 

The majority of the remnants are over luminous when 
compared to the observed Lx-Tx relation, while the 
most massive remnants are underluminous. The only 
remnants that appear to match the observed data well 



are the remnants of mergers between disks with 80% ini- 
tial gas fractions. 

The discrepancy for the most massive remnants ap- 
pears to be that they are simply underluminous, as this 
is the case for Lx-Tx as well as Lx-^sHowever, most 
of the other merger remnants appear consistent with the 
observed Lx~Lg relation, and thus it is unclear whether 
they are really overluminous or whether their tempera- 
ture is underestimated. The latter scenario could arise 
because of the radial temperature gradient present in 
most of the simulations. Hence, using a smaller aperture 
to measure the X-ray temperature or only detecting the 
central emission will yield higher Txthan that measured 
here. 

5.3. The a-Tx relation 

As a final comparison, we plot the X-ray tempera- 
ture against the stellar central velocity dispersion in Fig- 
ure ED As discussed in § 15.21 Tx is the mean X-ray 
weighted gas temperature and, as shown in ? 13.41 is a de- 
creasing function of radius. The central stellar velocity 
dispersion is measured within the half-mass radius from 
100 random projections. The typical la error owing to 
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Fig. 13. — Merger remnants from our entire series of gas fraction 
and size/mass simulations plot ted in the Lx-Lb plane. These are 
the same systems as in Fig. 1121 The black points arc galaxies from 
the X-ray catalog of lO'Sulliva n et alj I2q01al OFP01), and the 
crosses are X-ray luminous ellipticals from O'Sullivan ct al. (2003, 
OPC03). The dotted line is the estimated po int sou rce contribution 
log ^discrete = 29.5 log Lb from lCiotti ct al. (1991), and the dashed 
line is the best-fit relation for luminous ellipticals found by OPC03. 
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Fig. 14. — Bolometric X-ray luminosity vs. X-ray- weighted 
average temperature for our merger remnants. The crosses are X- 
ray-luminous ellipticals from OPC03, along with their best-fitting 
line. 



projection effects is ~10%. 

Provided that both the gas and stars are relaxed and 
in equilibrium, both a and Tx are a measure of the sys- 
tem potential. Also shown in Figure IT51 is a dashed line 
expected when the gas and stars are in equipartition, i.e., 



= kT, 



(4) 
is the 



where \i is the mean mass per particle, and m p 
proton mass. 

Figure 1151 plots a versus Tx for our series of merger 
remnants. In general, there exists a tight relation be- 
tween a and Tx for high mass remnants, signifying that 



Fig. 15. — Merger remnants and the central stellar velocity 
dispersion a vs. Tx plane. The crosses are X-ray-luminous ellipti- 
cals from OPC03. The dashed line marks equipartition of energy 
between X-ray gas and stellar velocity dispersion. 



these systems contain gas that is in hydrostatic equilib- 
rium throughout the remnant halo. There exists a trend 
for low-mass systems to peel away from the o~-Tx rela- 
tion. In particular, the remnants produced by gas-rich 
mergers show the largest discrepancy. These remnants 
are also strongly affected by the BH feedback (see, e.g., 
the amount of unbound gas mass in Fig. and their 
high temperature reflects this additional energy input. 

5.4. Building an X-ray halo 

Overall, our model for the production of X-ray emis- 
sion during the merger of two equal mass disk galax- 
ies successfully matches the rise and f all of X-ray emis- 
sion o bserved for the Toomre sequence l|Read &: Ponmar] 
1998). In addition, we find that the X-ray luminos- 
ity of the merger remnants span a wide range that de- 
pends on the initial disk gas fraction and size/mass. 
Compared to observed early-type galaxies of equivalent 
Lb, our models are broadly consistent with the data for 
Log(Ls)< 11 Lq, and underluminous for more massive 
remnants. However, because the hot gas is generally 
consistent with the Lx~Tx and u-Tx relations, we ar- 
gue that a disk-galaxy merger that includes the feedback 
from accreting supermassive BHs is a sufficient mech- 
anism to produce hot gas in the correct "state," i.e., 
metal-rich and in hydrostatic equilibrium, but in some 
cases there is simply not enough of this gas. 

As stated in the introduction, three mechanisms have 
been proposed to explain the growth of a luminous X- 
ray halo for elliptical galaxies, infalling tidal material, 
the reacquisition of material ejected in a galactic wind 
during the formation of the elliptical, and stellar mass 
loss. We discuss each of these in turn. 

5.4.1. Infalling tidal material 

One definitive conclusion from our simulations is that 
infalling tidal material contributes very little to the rem- 
nant X-ray emission. While tidal material continues to 
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rain down on the remnant long (~ 1 Gyr) after the nu- 
clei have coalesced, the mass of these cold tidal tails is 
small (~ 20%) compared to the mass of hot diffuse gas 
already in place. Not only is the mass in the tidal ma- 
terial not sufficient to significantly alter the X-ray lu- 
minosity, but the efficiency at which this infalling cold 
gas is converted to X-ray-emitting hot gas is less than 
30%. Thus, the diffuse hot gas is largely in place soon 
after the merger, as demonstrated in the T = 1.85 Gyr 
panel of Figure ^ an( i the tidal material serves only to 
"top off" the amount of hot gas. This is consistent with 
the theory that sh ells around ellipticals form from ac- 
creted tidal debris l)Hernauist & Spergellll992l) . because 
the stellar luminosity of these features is small compared 
with that of the host galaxies. 

5.4.2. Re- acquisition of ejected hot gas 

Within the context of the simulations we present, only 
the mergers that include BHs and produce a wind have 
a reservoir of ejected gas to reacquire. This is evidenced 
by Figures |3 and ^1 which show that the galactic wind is 
ubiquitous in small-mass halos and depends on the initial 
disk gas fraction. In some cases, the wind contains more 
mass than that which remains in the X-ray halo. 

Of course, all of our remnants will also accrete pristine 
gas through the hierarchical growth of structure, which 
we do not follow here. This gas is likely to be acquired by 
infall from the ambient intergalactic medium (which may 
be pre-enriched at a low level) or via discrete subunits. 
However, the high mctallicity observed in the halos of hot 
gas around local ellipticals ijHumphrev k, Buotel l2006() 
suggests that this is not the dominant mechanism to form 
the diffuse X-ray halo. Re-acquiring metal-enriched wind 
material allows for the growth of the X-ray halo while 
also maintaining a high metallicity. It is also possible 
that the enriched wind material will mix with the ambi- 
ent intergalactic medium, thereby diluting its metallicity. 
Effects such as this may be responsible for the scatter 
in the observed metallicities of X-ray halos or in radial 
metallicity gradients. 

5.4.3. Stellar Mass Loss 

Owing to the large population of stars formed dur- 
ing the merger, we expect a significant contribution to 
the hot gas content from ongoing stellar mass loss. The 
generation of diffuse hot gas from stellar mass loss has 
been discussed e xtensively in the literature (see, e.g., 
iCiotti et alJll991j) and is generally punctuated by a brief 
period of outflows triggered by ongoing Type I super- 
novae and stellar mass loss, followed by an extended pe- 
riod (~ 10 Gyr) in which gas is retained, and the hot 
gaseous halo then builds up. We emphasize that this 
"wind" is in addition to and would follow the galactic 
wind induced by BH feedback. This long timescale for 
the build up of the X-ray halo is consistent with the ob- 
servations o f elliptical galaxie s and m erger remnants, as 
discussed bv lO'Sullivan et all lj2001a|) . While our merger 
simulations do not alter this picture, they do provide a 
physically motivated initial condition that includes the 
mass of young stars. 

For the simulations we study here, the majority of the 
gas initially present in the progenitor galaxies is con- 
sumed by star formation. Thus the mass of young stars 



is typically 3-4 times that of the hot gas, which is ei- 
ther in the X-ray halo or unbound. While the amount 
of mass that is recycled from new stars back into the in- 
terstellar medium depends on the assumed initial mass 
function; if, for example, we assume a value of ~ 25%, 
then stellar mass loss will contribute an amount of hot 
gas that is roughly equivalent to what already exists from 
the merger alone. 

5.5. Galactic Winds 

A generic feature of the merger simulations that in- 
clude feedback from an accreting black hole is the pro- 
duction of a large-scale outflow of metal-enriched gas. 
This outflow carries a significant fraction of the avail- 
able gas mass and eventually deposits this gas through- 
out the halo and in some cases out of the galactic po- 
tential. Winds such as these may play a significant role 
in redistributing metals and may be required in order 
to terminate star formation and produce r ed ellipticals 
l|Springel et al.ll2005at IBopkins et alJl2006dl . 

Observationally, there is ev idence that winds su ch as 
these exist (for a review see iVeilleux et all 12001 . At 
low redshift, winds are mainly observed in star burst 
galaxies via metal-line absorption and emission stud- 
ies. It has been shown that there is a correlation be- 
tween the outflow rate and the star formation rate and 
that some of these winds are powerfu l enough to es- 
cape the galaxy l)Rupke et all l2005altJ: iMartirJ 12005^ . 
At high redshift, evidence for large-scale outflows comes 
from blues hifted absorpt i on lines, redshifted Lyq emis- 
sion lines llPettini et all l2002t lAdelberger et alJ l2003t 
Shapl ev et al.ll2003|) . and the excess of CIV absorption 
systems and a deficit of neutral hydr ogen near Lyman 
break galaxies (|Adelberger et al]l2005T) . 

It is currently unclear how closely the winds produced 
in the simulations resemble observed winds. We note 
that understanding these winds in more detail may re- 
quire embedding our simulations in a diffuse intergalac- 
tic medium, as many of the observations are actually a 
signature of the interaction between the wind and the 
local diffuse medium. This fact should be kept in mind 
when interpreting the fractional wind material presented 
in Figure^] Theoretical models predict that all galaxy 
halos, including those of spirals, should contain some hot 
gas component, yet it is currently unclear if these hot ga s 
halos are observed l|Toft et alJl20 02; Bens on et alJ l20001. 
In any case, including diffuse gas will provide a backstop 
against which the outflowing gas will interact, causing it 
to thermalize some of its energy, slow down, and possi- 
bly decrease the amount of outflowing material. We have 
performed one simulation in which both initial disks in- 
cluded a diffuse gaseous halo with one-quarter the bary- 
onic mass. In this experiment the wind punched through 
the diffuse gas, and this resulted in very little difference 
from the case without any diffuse gas. More work will 
be necessary to determine how the wind interacts with 
diffuse gas and to compare these simulations to the ob- 
servations. 

As a final point, we note that the absence of these out- 
flows in simulations without black holes does not mean 
that winds are only produced by accreting black hole 
feedback. This is quite apparent observationally, where 
systems such as M82 contain significant outflows and a 
corresponding star burst, but there is no evidence for an 
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accreting black hole. In our simulations the outflows are 
produced by point-like energy injection and one could 
easi ly consider other formulatio ns of star formation (see 
e.g.. lSpringel fc Hernauisti2003t) that would also produce 
galactic outflows. 

6. SUMMARY AND CONCLUSIONS 

Our hydrodynamic simulations show that collisions be- 
tween gas-rich disk galaxies generate appreciable X-ray 
emission. This emission comes from hot diffuse gas that 
is produced by strong shocks attending the merger pro- 
cess. The presence of an accreting black hole increases 
the X-ray luminosity during the merger through the pro- 
duction of a hot, metal-rich galactic wind. This galactic 
wind transports a significant fraction of metals from the 
galactic center to the diffuse gaseous halo and outside 
the potential well of the remnant galaxy. 

The wind-induced metal enrichment appears to be nec- 
essary to reprod uce the near solar met allicity of ob- 
served ellipticals (Hu mphrey fc Buotel 20061. Transport- 
ing metals out of the galaxy potential well may also help 
explain the high metallicity of the intra-group and intra- 
cluster media and the apparent deficit of baryons in lower 
mass galaxies. 

Our merger remnants span a wide range of luminosi- 
ties, yet are broadly consistent with X-ray properties of 
elliptical galaxies. This result suggests that galaxy merg- 
ers are a viable mechanism for generating the observed 
hot gas present in most interacting and merging systems. 
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In addition, we argue that both stellar mass loss and the 
reacquisition of galactic wind material are viable mech- 
anisms to increase the luminosity of the diffuse X-ray 
halo. 

This paper shows that galaxy mergers may play an 
integral role in the production of the hot, metal- rich 
gaseous halos that are present in most luminous ellip- 
tical galaxies. Using simulations similar to ours, recent 
studies demonstrate that gas-rich galaxy mergers that 
include black hole accretion and feedback can reproduce 
the black hole mass- ce ntral velocity disp ersion relation 
llDi Matteo et alJ 120051) an d the colors (|Springel et alJ 
2005a), fundamental plane ilRobertson et alJl2005|) . and 
kinematics of elliptical galaxies ijCox et alJl2005af . To- 
gether, these results pro vide strong evidence for a "cos- 
mic cycle" (see Fig. 1 of Hopkin s et al.ll2006al) of galaxy 
formation and evolution, in which gas-rich galaxy merg- 
ers and the growth of supermassive black holes play key 
roles. 

We thank Ewan O' Sullivan, Duncan Forbes, and Jes- 
per Sommer-Larsen for useful comments, and the ref- 
eree for suggestions that significantly improved this work. 
This work was supported in part by NSF grants ACI 
96-19019, AST 00-71019, AST 02-06299, and AST 03- 
07690, NASA ATP grants NAG5-12140, NAG5-13292, 
and NAG5-13381. The simulations were performed at 
the Center for Parallel Astrophysical Computing at the 
Harvard-Smithsonian Center for Astrophysics. 
\fCES 

Jenkins, L. P., Roberts, T. P., Ward, M. J., & Zezas, A. 2005, 

MNRAS, 357, 109 
Jog, C. J. & Solomon, P. M. 1992, ApJ, 387, 152 
Kim, D. & Fabbiano, G. 2003, ApJ, 586, 826 
Malin, D. F. & Carter, D. 1983, ApJ, 274, 534 
Martin, C. L. 2005, ApJ, 621, 227 

Mathews, W. G. & Brighenti, F. 2003, ARA&A, 41, 191 
McDowell, J. C, Clements, D. L., Lamb, S. A., Shaked, S., Hearn, 

N. C, Colina, L., Mundell, C, Borne, K., Baker, A. C, & 

Arribas, S. 2003, ApJ, 591, 154 
Mihos, J. C. & Hernquist, L. 1996, ApJ, 464, 641 
Mulchaey, J. S., Davis, D. S., Mushotzky, R. F., & Burstein, D. 

2003, ApJS, 145, 39 
Mushotzky, R. F., Loewenstein, M., Awaki, H., Makishima, K., 

Matsushita, K., & Matsumoto, H. 1994, ApJ, 436, L79 
Naab, T. & Burkert, A. 2003, ApJ, 597, 893 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1995, MNRAS, 
275, 720 

Nolan, L. A., Ponman, T. J., Read, A. M., & Schweizer, F. 2004, 

MNRAS, 353, 221 
O'Sullivan, E., Forbes, D. A., & Ponman, T. J. 2001a, MNRAS, 

328, 461 
— . 2001b, MNRAS, 324, 420 

O'Sullivan, E. & Ponman, T. J. 2004a, MNRAS, 354, 935 
— . 2004b, MNRAS, 349, 535 

O'Sullivan, E., Ponman, T. J., & Collins, R. S. 2003, MNRAS, 340, 
1375 

Pettini, M., Rix, S. A., Steidel, C. C, Adelberger, K. L., Hunt, 

M. P., & Shapley, A. E. 2002, ApJ, 569, 742 
Quinn, P. J. 1984, ApJ, 279, 596 
Raymond, J. C. & Smith, B. W. 1977, ApJS, 35, 419 
Read, A. M. 2003, MNRAS, 342, 715 
— . 2005, MNRAS, 359, 455 

Read, A. M. & Ponman, T. J. 1998, MNRAS, 297, 143 
Robertson, B., Cox, T. J., Hernquist, L., Franx, M., Di Matteo, T., 

Hopkins, P. F., Martini, P., & Springel, V. 2005, ApJ, 641, 21 
Robertson, B., Yoshida, N., Springel, V., & Hernquist, L. 2004, 

ApJ, 606, 32 

Rupke, D. S., Veilleux, S., & Sanders, D. B. 2002, ApJ, 570, 588 
— . 2005a, ApJ, 632, 751 



16 



Cox et al. 



— . 2005b, ApJS, 160, 115 

Sanderson, A. J. R., Ponman, T. J., Finoguenov, A., Lloyd-Davics, 
E. J., k, Markevitch, M. 2003, MNRAS, 340, 989 

Schweizer, F. 1980, ApJ, 237, 303 

Schweizer, F. & Seitzer, P. 1992, AJ, 104, 1039 

Schweizer, F., Seitzer, P., Faber, S. M., Burstein, D., Dalle Ore, 
C. M., & Gonzalez, J. J. 1990, ApJ, 364, L33 

Shapley, A. E., Steidel, C. C, Pettini, M., & Adclbergcr, K. L. 
2003, ApJ, 588, 65 

Sivakoff, G. R., Sarazin, C. L., & Carlin, J. L. 2004, ApJ, 617, 262 

Smith, B. J., Struck, C, & Nowak, M. A. 2005, AJ, 129, 1350 

Springel, V. 2000, MNRAS, 312, 859 

— . 2005, MNRAS, 364, 1105 

Springel, V., Di Matteo, T., & Hernquist, L. 2005a, ApJ, 620, L79 
— . 2005b, MNRAS, 361, 776 

Springel, V. & Hernquist, L. 2002, MNRAS, 333, 649 



— . 2003, MNRAS, 339, 289 

Strickland, D. K., Heckman, T. M., Weaver, K. A., Hoopes, C. G., 

& Dahlem, M. 2002, ApJ, 568, 689 
Toft, S., Rasmussen, J., Sommer-Larsen, J., & Pcdcrscn, K. 2002, 

MNRAS, 335, 799 
Toomre, A. 1977, in Evolution of Galaxies and Stellar Populations, 

p.401 

Toomre, A. & Toomre, J. 1972, ApJ, 178, 623 

Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARA&A, 43, 
769 

Voit, G. M. 2005, Rev. Mod. Physics, 77, 207 
Weil, M. L. & Hernquist, L. 1993, ApJ, 405, 142 
Xia, X. Y., Xue, S. J., Mao, S., Boiler, T., Deng, Z. G., & Wu, H. 
2002, ApJ, 564, 196 



